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The transport properties of finite nanotubes placed in a magnetic field parallel to their axes are 
investigated. Upon including spin-orbit coupling and curvature effects, two main phenomena are 
analyzed which crucially depend on the tube's chirality: i) Finite carbon nanotubes in a parallel 
magnetic field may present a suppression of current due to the localization at the edges of otherwise 
conducting states. This phenomenon occurs due to the magnetic-field dependent open boundary 
conditions obeyed by the carbon nanotube's wave functions. The transport is fully suppressed above 
threshold values of the magnetic field which depend on the nanotube chirality, length and on the 
spin-orbit coupling, ii) Reversible spin polarized currents can be obtained upon tuning the magnetic 
field, exploiting the curvature-induced spin-orbit splitting. 
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I. INTRODUCTION 

Carbon nanotubes (CNTs) have been the focus of in- 
tense study in the last decades as they exhibit remarkable 
properties that make them good candidates for molecular 
electronic devices. ^'^ Lately their application as build- 
ing blocks in spintronics has also been addressed. 
Nonetheless the importance of spin-orbit coupling in 
these systems has been for a long time underestimated. 
The spin-orbit interaction (SOI) was believed to be very 
weak, due to the low atomic number Z = 6 of car- 
bon. As a result spin degenerate states were assumed. 
Only recently it was proven experimentally^ that in the 
spectrum of ultraclean CNTs the effects of this coupling 
between spin and orbital degrees of freedom are clearly 
visible. This observation is in agreement with previous 
theoretical predictions^^^^ which argued that SOI could 
be significant in CNTs owing to their curvature. Un- 
derstanding the effects of this coupling on the spectrum 
of CNTs is essential for the successful manipulation of 
the different degrees of freedom of these systems. Sig- 
natures of spin-orbit coupling in the magnetic field de- 
pendence of the spectrum have been further identified 
in measurements^^ of relaxation and dephasing times in 
a CNT quantum dot and more recently in magnetocon- 
ductance experiments^^ performed on open CNT-wires, 
showing a double-peak feature caused by a spin splitting 
of the conduction bands due to the SOI. The observed 
spin-orbit energy splitting is compatible with previous 
experimental and theoretical values^' ^^'^^'^^ and further 
supports the relevance of spin-orbit effects in these sys- 
tems. Moreover new transport measurements^^ in CNT 
quantum dots with magnetic fields both parallel and per- 
pendicular to the nanotube serve to analyze the influence 
of SOI also in the cotunneling regime. 
The application of a parallel magnetic field modulates the 
distribution of the electronic spectrum of nanotubes, ^^'^^ 
as has been confirmed experimentally in several opti- 
cal and transport measurements. ^^^^ Further theoreti- 
cal studies^^^^ reveal that the transport properties are 



directly related to the magnetic-field modulated density 
of states. The SOI -modified band structure and its im- 
plications on the electronic properties of CNTs in a mag- 
netic field has been in the focus of extensive theoreti- 
cal work^^'^^^^. However, at present only few theoreti- 
cal works exist^^'^^ which investigate quantum transport 
through CNTs including magnetic field, curvature and 
SOI effects. 

In this paper we analyze the transport properties of 
finite-size carbon nanotubes in a parallel magnetic field 
in the regime of strong nanotube-lead coupling, using a 
tight-binding Hamiltonian for the nanotube system and 
Green function techniques. A magnetic field parallel to 
the tube axis is included via the Aharonov-Bohm phase 
and Zeeman effect, and we take into account the spin- 
orbit coupling, which is responsible for the appearance 
of a spin-splitting and of spin-flip processes. The spin- 
splitting may allow for preparation of states with con- 
trolled spin and valley (the so-called isospin) degrees of 
freedom. The inclusion of the SOI and the magnetic 
field leads furthermore to the breaking of the electron- 
hole symmetry, as we notice in our analysis of these sys- 
tems. Without an electron-hole symmetry spin-up and 
spin-down states may cross the Fermi energy at differ- 
ent magnetic fluxes giving rise to peaks in the magne- 
toconductance corresponding to spin-polarized currents 
around the Fermi energy similar to those experimentally 
observed in Ref. 14 and as identified in the present work. 
For zigzag and chiral nanotubes the finite size of the sys- 
tems leads to a localization of the states in the magnetic 
field with a consequent suppression of current. ^^'^^ We 
show that, as a consequence of the SOI , this localization 
starts at a threshold value of the magnetic field which is 
different for the two spin species. This value can be fur- 
ther tuned with the chirality, radius or length of the tube, 
allowing for controlled polarized states. Furthermore we 
find that numerical transport calculations in finite-size 
CNTs are to a very good extent reproduced by using 
an analytical model Hamiltonian in the reciprocal space. 
The proper boundary conditions are crucial for the un- 
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derstanding of finite-size effects. Tiie analytical approach 
provides as well a fundamental insight into the nature of 
the states taking part in the transport processes. 
Part of the discussion on localization was presented in 
abridged form in Ref. 38. The scope of the present work is 
to analyze various consequences of SOI on the transport 
characteristics of CNTs strongly coupled to leads. 
The paper is organized as follows. In Sec. II we introduce 
the tight-binding approach that we use for modelling our 
system, where we keep the interactions to nearest neigh- 
bors only because, despite its simplicity, this approxima- 
tion is known to provide a good agreement with experi- 
ments.^ Our original tight-binding Hamiltonian is mod- 
ified to consecutively include the effects of an applied 
axial magnetic field, the curvature of the nanotube and 
the spin-orbit coupling. Starting from this generalized 
tight-binding model for the real-space description of the 
CNTs, an effective Hamiltonian for the tt bands is ob- 
tained, which is formally the same as the one derived pre- 
viously by other groups. ^^'^^ This Hamiltonian will allow 
us in Sec. HI to make predictions for the energy spectrum 
when taking into account the boundary conditions to be 
satisfied at the CNT ends. Finally the transport calcula- 
tions are shown and discussed in Sec. IV- VI. An analysis 
of the contribution to the CNT conductance properties 
of the different effects under consideration is carried out 
in Sec. IV, followed by a comparison of numerical and an- 
alytical results and an examination of the localization in 
the magnetic field of previously extended states in Sec. V. 
The great potential of these CNT systems for their use 
in spintronics is demonstrated in Sec. VI, where the spin- 
polarized currents at small bias voltages are investigated. 
A summary of the main results is found in Sec. VII. 

II. MAGNETIC FIELD AND SPIN-ORBIT 
COUPLING EFFECTS 

We wish to describe a finite single- wall carbon nanotube 
with a chiral vector^^ C as sketched in Fig. 1. The chiral 
vector uniquely defines the nanotube and determines its 
unit cell, radius R and other properties such as its chiral 
angle 7^, which is defined such that for zigzag CNTs 77 = 0. 
Neglecting for the moment curvature effects and spin- 
orbit interaction, we restrict ourselves to one 2p orbital 
(pz) per site, as these orbitals give rise to the tt molecular 
orbitals, which in these systems yield the major contri- 
bution to conduction properties. Then the Hamiltonian 
describing our system can be written in the tight-binding 
approach as follows: 

Ho = ^ e2p4c. + t%c\c., (1) 

i (id) 

where the indices (z, j) indicate nearest neighbor atom 
sites and the summation is extended over all the points 
in the lattice. The onsite energies are e2p and t^j are the 
hopping parameters between neighboring sites. We will 
shift our energy scale in order to have vanishing on-site 



energies, setting = 0. The parameter t^j is given by 
the hopping between nearest neighbor pz orbitals, Vpp. 
The above Hamiltonian in real space needs to be modi- 
fied in order to incorporate the effects we are interested 
in: (A) a parallel magnetic field, (B) the curvature of 
the carbon nanotubes, (C) the spin-orbit interaction. In 
the following we will therefore find the modified hopping 
parameters providing these effects. 

From now on we shall be using two sets of coordinates 
to describe our system. The first, (x^,?/^,^^), denotes 
the local direction of the orbitals at the atom i of the 
nanotube. The direction Xi is tangent to the circumfer- 
ence of the tube, yi parallel to the nanotube axis, and Zi 
is always perpendicular to the nanotube surface as seen 
in Fig. la. The second system, (X, F, Z) accounts for 
the global cartesian coordinates of the nanotube, with 
the y-axis coinciding with the tube axis. In the global 
cylindrical coordinates (i?, 6>, Y) the angle 6 is measured 
from the positive Z-axis and R is the nanotube radius. 
Therefore when we consider the atomic orbitals, the tt or- 
bital is in the Zi direction, perpendicular to the nanotube 
surface and the magnetic field will be applied along the 
y-direction. 



A. Parallel magnetic field 

In the presence of a magnetic field the Hamiltonian 
describing our system is obtained by the substitution 
— iW — iW — f A, where A is the vector potential 
of the magnetic field B. As shown by Luttinger,^^ the 
effect of the magnetic field in the Schrodinger equation 
is greatly simplified if the Bloch functions are modified 

by a phase factor exp (if /^^ A(rOdr') .^^ Therefore each 

matrix element in the Hamiltonian is obtained by mul- 
tiplication of the corresponding matrix element in zero 
magnetic field by a phase factor: 

4(B)=4(0)e"^^^'^^^^ (2) 

known as Peierls phase factor, where the integral is 
taken along a line joining the sites i and j. 
In the case of a uniform static magnetic field in the di- 
rection of the nanotube axis, B = (0,50,0), the vec- 
tor potential can be given by A = ^(— Z, 0,X), or 
(0, ^Bor, 0) in cylindrical coordinates, where we have 
taken the Coulomb or transverse gauge. With this choice 
of gauge the Peierls phase factor is given by 

e^fe;Bor'>ie ^^i^(9.-9,)^ (3) 

where (j) is the magnetic flux in the cross section of the 
nanotube perpendicular to its axis, = h/e is the flux 
quantum, and Oi — Oj is the angular difference between 
the azimuthal coordinates of sites i and j. 
The Zeeman effect appearing in the presence of a mag- 
netic field is taken into account by including in our Hamil- 
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FIG. 1: System and coordinates, a) Schematic illustration of 
the carbon nanotube with the local and global coordinate sys- 
tems used for its description. The direction of the magnetic 
field is parallel to the tube axis, b) Projection of the nan- 
otube onto the X — Z plane, i and j schematically show the 
projection of two nearest neighbor atoms (at an exaggerated 
distance in order to visualize the different variables). The an- 
gles are measured from the positive Z axis and \0j — | is very 
small for large radii, so that in fact both angles can be seen as 
0. R is the radius of the nanotube, whereas Hij is the vector 
connecting the two neighboring atoms (here we plot actually 
the projection of this vector), c) Real space representation 
of a fragment of a graphene lattice with translational vectors 
ai and a2. The chiral vector is C and the unit vectors x± 
and X|| reflect the parallel and perpendicular directions with 
respect to the nanotube axis. The three B nearest neighbor 
atoms of an atom A are located at positions given by with 
/ = 1, 2, 3. d) Reciprocal space for the graphene, spanned by 
the vectors bi and b2. The central hexagon is the symmetric 
Brillouin zone. Its corners are the so-called K and points 
or Dirac points. 



tonian the term Hz = ^M^B • a, where ge is the elec- 
tronic g-factor, /iB = ^ is the Bohr magneton and a 
is Pauli's spin matrix vector. With the magnetic field 
directed along the F-axis, only the Pauli matrix ay re- 
mains, which conserves the spin in the F-direction with 
eigenvalues ±1 for parallel/antiparallel spins in this di- 
rection. Expressing the magnetic field Bq in terms of 
the magnetic flux and with the g-factor set to 2, the 
Zeeman term is then 



nz = 



6o mK^ 



cry, 



(4) 



giving a corresponding shift of for spins 

up/down in the y-direction. 



Due to the curvature of the nanotube lattice, the 2pz or- 
bitals which are parallel to each other in the graphene 
are not parallel anymore (see Fig. 2). This new inter- 
play of the chemical bonds has to be taken into account, 
and we have to recalculate the hopping between pz or- 
bit als accordingly. This hopping will now have not only 
components due to orbitals normal to the atomic plane, 
but also due to tangent ones. The hopping between two 
neighboring pz orbitals at sites i and j is then 

{Zi\n\zj ) = V;^ n^{zi)-n^{zj) + V;^ nll(^0-n"(^i) (5) 

where 11^(2:^) and nll(z^) are the normal and tangential 
components of a unit vector in the direction of the Pz 
orbital at site i, ), and V^^ and V^p are the transfer in- 
tegrals giving rise to a and tt orbitals in flat 2D graphene, 
respectively. We use throughout the paper round brack- 
ets to indicate states in the absence of spin-orbit inter- 
action and we will leave the normal bracket notation for 
the states that include this perturbation, as introduced 
in Ref. 10. The parallel component of n(a^) is its pro- 
jection along Hji given by nll(ai) = ^^^^^^5^ Rj^, and 

I i I 

n^(a^) = n{ai) — nll(a^), where Hji = Hj — is the 
vector connecting the two sites. Therefore 

{zi\n\zj) = V^p n{zi)-n{zj) 



For the case of Pz orbitals this means that 

i zi\n\z, ) = v;^ cos i0i-0,) 

-{v;^-v;^)-^[i-cos{e,-e,)]\ 

where ac = \Rji\ is the distance between neighboring 
atoms. 
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C. Spin-orbit coupling (SOI) 

Carbon 2p orbitals present a weak spin-orbit interaction 
(SOI ).^^ Because the effect of this interaction is small we 
shall treat it as a perturbation, and consider its effects 
to the lowest order. 

We follow here the method used by Ando,^^ and adapt it 
to the case of spins in the direction of the parallel mag- 
netic field. First we have to calculate the modification 
of the 2pz atomic orbitals once the SOI is included. The 
Hamiltonian describing this interaction is given at the 
atomic level by^^ 

^so = ^{^Vxp)- a, (8) 

I 



where V{t) is the atomic potential, m is the free-electron 
mass and p is the momentum operator 

As described in Appendix A we include the SOI Hamil- 
tonian as a perturbation and obtain to first order the 
modified Pz orbitals, \zjs) with s =ty^iy^ denoting the 
spin up/down in the F-direction. The perturbed orbital 
\zjs) does not have anymore a well defined spin and the 
spin in the axial direction s indicates merely the main 
character of the state. Moreover \zjs) contains also an 
admixture of the p orbitals building the a bonds. 

The transfer integrals for these perturbed states are now 
given by 



= {z, ty I n \z, t,) ^ v;^ cos [Oi - ej) - {v;^ - v;^)-^ [cos {0, -0,)- if 

+ 2iS sin {0, -0,) + {V;^ - V7p) J sin {0^ - 0j) [1 - cos {0, - 0,)] | , 

hi = {Zi iy I n \Zj iy) = {Zi I H \Zj ty)*, 

tn = % I -H \z, iy) -5 {e-''^ + e-'^O (T/^ - V;^)-^ [cos {0, - 0,) - 1] , 

ht = i^i iy I ^ iy) = -(^i iy I ^ V3 iyT ^ 



where Yji is the distance between the interacting atoms 
projected in the axial direction and (5 is a dimensionless 
parameter measuring the SOI strength as defined in Ap- 
pendix A. 

With these recalculated hopping parameters we can build 
our Hamiltonian in real space and numerically compute 
the transport properties of our systems, c/. Sec. IV. In 
order to gain a deeper insight in the results, it is im- 
portant though to approach the problem analytically as 
well. Hence in the coming section we derive from Eq. (1) 
and with Eq. (9) a low energy Hamiltonian in k-space. 
In this way we will not only be able to understand the 
observed numerical features but also to make analytical 
predictions and save computational costs. 



D. Low energy Hamiltonian 



evaluated in Appendix B explicitly in the limit of large 
radii, so that ac/R is a small quantity, and in the limit 
of low energies, such that an expansion of the Hamilto- 
nian around the K and points of the reciprocal lattice 
(see Fig. Id) is justified. The matrix block H^^(k) be- 
comes in this latter limit H^^(k), valid for small vectors 
K around the Dirac point rK with r = +1 for K and 
r = -lfor (k = rK + A^). 

With the definitions A/c^ = f^rp cos 3r], Aku = 



-^p' sin 37?, Akl^ 
with p = 1 



25p 
R ' 



and A 



flip 



pp PI 
4i? VJ„ 



8 K 



pp pp ry^l 



^YIv_X^^ we can write 



Wab{i^) = ^t^Fe"^^^J(r/^± +rAfcl) -i(/^|| +rAfcf| 



The Hamiltonian of carbon nanotubes can be written 
in the basis of the Bloch functions for the A and B 
sublattices in the same way as for graphene:^^ ^(k) = 

...i^^'^'f^) ^AB(k) = Ef=ii<e>k-d<. As 

shown in Fig. Ic, the vectors connecting neighboring 
atoms are denoted by d^. As the spin degeneracy is bro- 
ken, the transfer integral ti in (k) becomes now a 2x2 



matrix 



These matrix elements are 



hvpe 



where vf 



rAk 



so 



-ay — hvpe 



(10) 



h ' 



8.6 • lO^m/s is the group velocity 
at the Fermi points. Notice that the spin-flipping term 
is actually a rotation into the local coordinate system 

(-e^^~^o^) i (cos 6ax -sinOaz^ = icr^,,, ^ij being 
the direction tangential to the tube at the middle point 
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of the bond between i and j atoms. 
The eigenvalues of are then 



E± (r, a, /^x, ^ii) = ±hvF {{ri^i_ + rAklf + {Akf^)^ 



+ 2a [((A^!^)2 + (a^ii + rA^P") A 
+ (r,.x+rAfci)^(Afei^n'^'}'^', 



2 

/Hp 



(11) 



with a = ±1 referring to the spin, so that we could merge 
all four eigenvalues into this expression. 
If we neglect the term allowing for spin flip (AfUp)^ which 
is small in comparison to other contributions, we have 



E± (r,a,A^Y,/^j|) =±n^;Fy(A^Y)'+ , (12) 

equivalent to the dispersion relation for the planar 
graphene, but with the definitions 



k'w = K\\ + TAk?\. 



(13) 



Neglecting the spin-fiipping terms we recover a well- 
defined spin in our problem, so that we can distinguish 
between spin- up and spin-down eigenfunctions. 
The magnetic field modifies the dispersion relation fur- 
ther via the Aharonov-Bohm term, which changes f<i± by 



so that K.^j_ in Eq. (13) becomes 



/^x + Akl + crA/ci^ + T- 



(14) 



The contribution of the spin-orbit interaction on the di- 
agonal of the matrix blocks can be also interpreted as an 
effective magnetic field and we can write Ak^^ = 
with ^50/^0 = 2(^p. The other contributions, Akj_ and 
A/cy, arise from the consideration of the curvature in the 
structure of the nanotubes, as it can be seen in Ref. 38. 
The Zeeman effect modifies further the dispersion rela- 
tion given in Eq. (12) by introducing a spin-dependent 
shift: 



E± (^r,cr, = ±hvF\l {i^'^f + ^cr/J^B |B| , 

(15) 

that is, levels of up (down) spins are correspondingly 
shifted up (down) in energy. Ref. 30 introduces an ad- 
ditional effective valley- dependent Zeeman term which 
leads to a small energy splitting in the Dirac cones and 
introduces a further electron-hole asymmetry. The mag- 
nitude of this shift only slightly modifies this spectrum, 
and will therefore not be considered in the following cal- 
culations. 



III. FINITE SIZE NANOTUBES 

Eq. (15) allows us to calculate analytically the spectrum 
of CNTs in a parallel magnetic field, including SOI and 
the Zeeman effect, when in addition proper boundary 
conditions are imposed that yield the appropriate quan- 
tization of the vector k,. 

First, as a CNT is obtained by rolling a graphene layer 
into a tube, in the angular direction its wave functions 
always obey periodic boundary conditions: (C) = 
exp (ik • C) tjj (0) = e^'^^'^ip (0). This boundary condition 
leads to a quantization of the transverse component of 
k, k± = with l± being a positive integer. We 

can obtain the energy spectrum of infinitely long carbon 
nanotubes from the dispersion relation of graphene by 
restricting the allowed values of k by the quantization 
condition given above. The quantized k± are then ob- 
tained from the relation = k±— rK^, which together 
with Eq.(14) determines the quantized values. 
To obtain the quantization of /^y, and therefore of for 
CNTs of finite length we have to impose open boundary 
conditions in the axial direction, that is, we impose that 
our wave function vanishes at the ends of the tube. To 
this extent we have to consider the underlying structure 
of the nanotube. Here we distinguish between two cases: 
i) armchair nanotubes, where each end contains the same 
number of atoms from the A and B sublattices and ii) 
zigzag and chiral nanotubes, where each end is formed 
predominantly by one type of atom. Because the same 
boundary conditions apply for zigzag nanotubes and any 
other chiral CNT,^^ provided the nanotube edge is a so- 
called minimal boundary (there are no atoms with only 
one neighbor), we shall derive them for the zigzag case. 
In zigzag CNTs one end is entirely composed of A sub- 
lattice atoms and the opposite end is formed only by B 
sublattice atoms. 

Zigzag nanotubes: The open boundary conditions in this 
case imply that the wave function must vanish on one 
end at the missing A atoms and on the other end at the 
missing B atoms. We therefore impose 



i^rA {{X±,X\\ = L), E±) = 0, 
i^rB ((^±,^11 = 0),E±) = 0, 



(16) 



where the wave functions V^rA, V^rB are the sublattice 
components of iIj^- (r, £^±), the eigenstate at a given en- 
ergy E{n'^^ n'^^) and at the Fermi point rK. Because the 
solutions with ti'^^ and AC 1 1 are degenerate in energy, as 
can be seen from the energy dispersion, Eq. (15), the 
eigenstate satisfying Eq. (16) will be a linear combina- 
tion of the eigenstate for and the one for — or 
equivalently the one for n\\ and — 2rA/cy . As shown 
in Ref. 38 the condition to be fulfilled by the k values is 
then 



A2K'n L 



(17) 
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This equation has both extended and locahzed solutions 
and a trivial solution for = 0. If is real the condi- 
tion given by Eq. (17) becomes 



TKjj^ — ^Cii COt^Cii-Zv, 



(18) 



which represents the extended solutions. Otherwise if k,'^^ 
is a pure imaginary number Eq. (17) becomes 



= Im(/^j|) coth (lm{K^^^)L^ . 



(19) 



As /^ii is imaginary, the solution to this equation describes 
evanescent waves. Eq. (19) has two non- trivial solutions: 
> 1/L for the K-point and k,'^ < —1/L for the K'- 
cone. 

For =0, Eq. (17) yields as allowed values of /^n: 



k|| = (2n+l)— , 



(20) 



with n G Z. 




imaginary, K' 
real, K' 
imaginary, K 
real, K 
trivial 



FIG. 3: Allowed Kj_ and Atji values for a zigzag- like CNT 
fulfilling the boundary condition given by Eq. (17). At = 
0, as expected, — + ^ with n G Z and when including 
the spin the Kramers degeneracy holds. 



The allowed solutions of Eq. (18) and Eq. (19) are plotted 
in Fig. 3. 

Armchair nanotuhes: The infinite armchair CNT will be 
always metallic in the 7r-orbital approximation, as its chi- 
ral angle of 30° assures the cut of the K and points. 
Bloch waves on branches of the linear dispersion relation 
with positive (negative) group velocities, vp^ are called 
right (left) movers, as seen in Fig. 4. 
The armchair case was not discussed in Ref. 38. Here 
we extend the treatment of Ref. 50 to include SOI and 
magnetic field effects. To impose the open boundary con- 
ditions we have to build once more linear combinations of 
Bloch waves (in the expansion around the Fermi points) 
which are degenerate in energy. In the presence of an off- 
set A, which accounts for the case that |K|| | = J (m — A) 
with m being a positive integer, the branches that are 




FIG. 4: Linear dispersion relation of finite size armchair 
carbon nanotubes with open boundary conditions, a) The 
band structure of an armchair CNT with periodic boundary 
conditions around the two inequivalent Fermi points, K and 
K^ For finite- length tubes only certain values of are al- 
lowed which are here indicated by the sketched vertical cuts. 
Branches of left movers are indicated by / whereas r indi- 
cate those of right movers, b) These four branches can be 
combined into the band structure shown on the right, which 
represents the solutions for the open boundary problem. We 
have chosen explicitly the case where z^:^ = and the allowed 
At 1 1 do not lie symmetrically around the Fermi points, in order 
to easily observe the bands of equal energy. Including SOI 
and the magnetic field in k'^ leads to shifted cuts of the Dirac 
cones and we would have the same picture as here but instead 
with parabolas. 



equivalent in energy are the r(/)-band at K and the l{r)- 
band at K^ The linear combinations of wavefunctions 
are then: 

i^r (r,A^) = ^ (<I>^'(r,(A^^,/^||)) -<I>f^(r,(A^^,-/^||))) , 



(r, 1^) 



^($f'(r,(/.^,/.||))-$^(r,( 



V2 



(21) 



where ^^j^ are the Bloch wave functions for the l/r 
branches. We have therefore to mix states which belong 
to different Dirac cones. 

An analysis of the wave functions ip± reveals that they 
are proportional to sin ((Ky + x\\). The open bound- 
ary conditions on these wave functions are given by^^ 



0, 
0, 



(22) 



and require the wave functions to vanish on atom 
sites belonging to both sublattices (as opposed to the 
case of zigzag nanotubes). They are then fulfilled if 
(K|| + /^ii) L = 7rn, with n G Z. That is. 



{n^ + A) , with G 



(23) 



A value of A < 1 different from or ±^ causes a mis- 
match eA = hvpj^ (2A — 1) between the f and / branches 
as seen in Fig. 4. One of the bands will therefore be en- 
ergetically favored if an electron is added to the tube. 



To include the effects of SOI and an applied magnetic 
field it suffices to note that the Bloch wave functions 
building up ?/^~^^~ for a fixed spin a are still degener- 
ate when considering Eq. (15). Indeed n'^^ includes a 
valley-dependent shift, so that the Bloch wave functions 
<I>^^(r, (/^x, /^ll), cr) and $^(r, (a^^, — A^y), <j) are charac- 
terized by = A^ii — A/cy and = — A^y + A/cy and are 
therefore degenerate in energy. 

The allowed values for in the case of armchair CNTs 
are all real and no imaginary solutions are possible. Thus 
all states building the energy spectrum of armchair tubes 
will correspond to extended states. 



IV. TRANSPORT PROPERTIES IN A 
PARALLEL MAGNETIC FIELD 

The transport through the systems we studied in the pre- 
vious sections is calculated by connecting the CNTs to 
metallic leads, as schematically shown in Fig. 5. Nan- 
otubes of various chiralities and lengths are considered. 
Their transport characteristics are obtained from the 
transmission values gained by applying the Landauer 
formula, using Green function techniques. In particu- 
lar, we determine the elastic linear response conductance 
via the Fisher-Lee formula for the quantum mechan- 
ical transmission^^: G = ^Tr {Fl 5 Fr ^"f" } , where 
Fl/r = i(5]L/R - El/r is the self energy as- 

sociated to the left or right lead respectively, and Q is 
the Green function of the central region dressed by the 
electrodes, Q = {E — H — El — 5]r)~^. For simulating 
bulk metal electrodes we will consider wide band leads, 
i.e. ^wb{E) ^ -ilm5](£; = Ef). 



DoS 



DoS 




FIG. 5: Schematical drawing of the set-up of our system, 
where a suspended CNT connects the source and drain leads, 
between which a bias voltage Vb is applied. 

We will now show in detail the spectrum and transport 
characteristics of a finite zigzag nanotube of chirality 
(9,0) with a length of hundred unit cells (about 42 nm). 
In order to see the influence in transport of the different 
effects considered, we will include them gradually. More- 
over we will consider a rather large value of 5 to observe 
clearly the effects of the spin-orbit interaction. All other 
parameters used for the following calculations are taken 
from Ref. 54. 

With the help of the effective Hamiltonian and its eigen- 
values given by Eq. (15) we can see that for the bare 




FIG. 6: Density of states of a (9,0) CNT of a length of about 
42 nm (100 unit cells) a) within a vr-orbital tight-binding ap- 
proximation, neglecting all other effects, b) with the inclusion 
of curvature effects as described in Sec. II B, c) and taking into 
account the Zeeman effect. Finally in d) spin-orbit coupling 
is additionally considered. The value of the SOI parameter S 
is set to 0.013 and Y,wb{E) ^ -i0.85eV. 



di/dv 



di/dv 




c) 






0.2| 


> 











-0.2 
C 



dl/dV 



FIG. 7: dl/dV characteristics in units of the conductance 
quantum Go, for a (9,0) CNT of a length of about 42 nm as in 
Fig. 6 with the respective effects being successively included. 



TT-orbital model (Fig. 6a and Fig. 7a), there is a four- 
fold degeneracy at zero magnetic field: a spin- and a 
valley-degeneracy. The eigenvalues at 5 = correspond 
to the same allowed value of (as we are dealing with 
a metallic CNT, specifically i<i± =0), and different al- 
lowed values of as determined by Eq. (17). The valley- 
degeneracy is broken with the onset of the parallel mag- 
netic field (Fig. 6a and Fig. 7a). As argued in Ref. 55 
and demonstrated analytically below, cf. Eq. (24), this 
is the expected behavior if one assumes that i) the shift 
in energy due to the interaction of the orbital magnetic 
moment and the magnetic field is AE = —/J'^rb ' 
that ii) the orbital magnetic moment is reversed when 
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moving from one Dirac cone or valley to the other. This 
implies that the states in the K Dirac cone will move to- 
wards the Fermi energy as their orbital magnetic moment 
is parallel to the magnetic field for conduction electrons 
and antiparallel for valence electrons. In the cone 
the behavior will be inverted with states moving apart 
from the Fermi energy. This is indeed the evolution of 
the spectrum seen at low magnetic field. As the magnetic 
field is increased the evolution of the states becomes com- 
parable in both valleys: all states move away from the 
Fermi energy, both in the valence and conduction band. 
We can understand this behavior in both regimes by look- 
ing at the limits of Eq. (17). For small magnetic fluxes 
the components of the vector k' can be written as 
and n'u ~ ^(2^ + 1) + A/^m with n G Z. 



In this limit Eq. (17) yields A/^m = 



. The 



R 7r(2n+l) 0o 

spectrum given by Eq. (12) reads then to first order in 



.2L 



(2n + 1) 



1-r 



2r (j) 
L 4 



i?7r2(2n + 1)2^0 



(24) 



that is, we have a linear evolution with the magnetic 
fiux with slopes of different signs for the two indepen- 
dent Dirac cones. In other words we obtain the effec- 
tive orbital moment at low fields being \/J^orb\ — 7r(2nTi) ' 
which depends on the discrete value of /^y, and decreases 
as we move away from the band edges (corresponding to 
=0). This is a signature of the finiteness of the tube, 
which could not be captured by the classical estimate of 
the orbital magnetic moment done in Ref. 55. 
The finite values of described by Eq.(20) assure the 
presence of a gap in the absence of a magnectic fiux, as 
observed in Fig. 6a and Fig. 7a. 

At higher fields the evolution of the resonances with the 
magnetic fiux is again linear. Now the allowed values of 
remain nearly constant approaching asymptotically 
wi^h n > G Z, and the dominant contri- 
bution to the evolution of the states is given merely by 
that is, evolves linearly with a slope of as given by 
Eq. (14). Moreover this slope is independent of the val- 
ley degree of freedom. At intermediate fields a crossover 
between the two different linear regimes appears. 
As seen in Fig. 6b and Fig. 7b, the inclusion of the cur- 
vature modifies the spectrum of the CNT widening the 
gap at zero magnetic fiux. This can be understood by 
observing that (/^Y)^ larger after the inclusion of Akj_ 
(see Eq. (13)) and at the same time a larger value of l^^^l 



implies a larger 
spin-degenerate. 



(see Fig. 3). The states are still 
We also observe that the states now 



cross as they evolve with the magnetic fiux. The cross- 
ings are possible because of the curvature-induced shift 
Ak^ of the perpendicular momentum, which is valley 
dependent. We have therefore two different evolutions 



for /^^ for the K and cones, which allow crossings of 
states of different valleys. 

In Fig. 6c and Fig. 7c we introduce the Zeeman effect in 
this picture. The spin-degeneracy breaks at finite fields, 
and an additional evolution, opposite for spins up and 
down, competes with the evolution due to the orbital 
magnetic moment. 

Finally, SOI introduces a splitting of the spin species al- 
ready at zero magnetic field (Fig. 6d and Fig. 7d). Now 
tz^^ is shifted towards negative values for spin up and to- 
wards positive values for spin down. A detailed analysis 
of the momentum shift caused by the different effects 

combined with the diagram of solutions for 

(Fig. 3), allows us to conclude that now the degeneracy 
at zero field is between a K ^ and a t subband and 
between a K t and a | subband. This degeneracy is 
broken at finite magnetic fields. Including the Zeeman 
effect introduces a different evolution with the magnetic 
field for the two spin species. This leads to the break- 
ing of the electron-hole symmetry in the spectrum of the 
CNT as it can be seen in Fig. 6d. In the corresponding 
dl/dV characteristics we recover a symmetric picture, 
as we are assuming a symmetric voltage drop across the 
junction. The chemical potential of the right and left 
lead is set at Ep ± with V5 being the bias voltage. 
In this case at zero temperature the Landauer formalism 
allows us to write the differential conductance as 



dl _ 2e^ T {Ef , 2 
dV ~ 



-T{E, 



2 > 



(25) 



h 2 
which is symmetric in the bias independent of the trans- 
mission T. 

The Zeeman splitting leads as well to the breaking of 
the quantum fiux periodicity in the magnetic field, which 
has been checked numerically (not visible in Fig. 6 and 
Fig. 7 where a smaller range of the magnetic field has 
been chosen for clarity). 



V. LOCALIZATION IN THE MAGNETIC FIELD 

In Fig. 6 and Fig. 7 we observe a remarkable feature. The 
density of states acquires a pronounced peak at zero en- 
ergy above a threshold value of the magnetic field. Once 
the Zeeman effect is considered, this peak is split into 
two, appearing now at finite energies. At the same time 
we observe that the differential conductance for these 
states is strongly suppressed. This feature refiects the 
onset of the localized solutions^^'^^ described by Eq. (19). 
States with imaginary k,'^^ occur in the vicinity of the 
Fermi point K at = 1/L, and for instead at 
K,'^ = —1/L. The threshold value in the magnetic field is 



therefore given by 



R 



L 



• /^x - AA:^ - aAk^ 



so 



(26) 



where we have made use of Eq. (14). That is, the onset 
for the imaginary solution in the K cone is around R/L 
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component (see panels c and d of Fig. 8) . For the corre- 
sponding valence eigenstate localizing as a spin down we 
obtain the same picture with just a different value for the 
onset of the localization {(j)ioc ^ 0.029(/)o for r = 1 and 
cr = —1), so that for the magnetic flux of ^ = 0.0500, the 
localization at the edges is stronger. 
Many nanotubes (most notably zigzag CNTs) present 
"native" localized states, present at their edges even in 
the absence of the magnetic field. ^^^^ They correspond 
to states from higher energy subbands, with i<i±_ = 
where hc'^ fulfills the localization condition naturally. The 
native localized states are not seen in these figures be- 
cause of the rather strong coupling to the leads, but are 
recovered if we reduce it. 



FIG. 8: Spatially resolved wave function (obtained numeri- 
cally) of a (9,0) CNT with 100 unit cells and SOI correspond- 
ing to the valence eigenstate getting localized as a spin-up 
state (see Fig. 6). In the top figures, dots represent the atomic 
positions and their color corresponds to the amplitude of the 
wave function at the site (going from zero to its maximal 
value). The two sublattices A and B are plotted separately. 
The corresponding magnetic fluxes are a) = 0.02 00 and 
b) = 0.0500. The decay of the wave function amplitude for 
these two flux values is shown in c) and d) respectively, where 
sublattice A is plotted with a dashed line and a solid line is 
used for the results corresponding to sublattice B. Red marks 
the spin-up component of these states while blue is used for 
spin-down, which is negligible here. 

with a shift towards larger fields depending on the chi- 
rality (being maximal for zigzag CNTs), and a further 
shift depending on the SOI parameter 5 towards larger 
fields for spin up and towards lower fields for spin down. 
Analogously for K^ 

These eigenstates gradually get localized at the edges of 
the CNT and we can see in Fig. 8 the onset of this local- 
ization for the spin-up state of the valence band which 
localizes above the Fermi energy. Fig. 8a corresponds 
to a magnetic flux below the threshold value (j)ioc and 
the magnetic flux in Fig. 8b is already above this value, 
which for this CNT is 0^oc ^ 0.04300 (for r = 1 and 
a = 1). We therefore see in the decay of the wave func- 
tion amplitude (Fig. 8c and Fig. 8d) its change towards 
a localization at the edges with an exponential decay for 
the larger magnetic field value as it corresponds to an 
imaginary momentum. The A and B components still 
overlap for this magnetic field, so that the corresponding 
conductance is reduced but not suppressed yet. In both 
cases they are eigenstates of the angular momentum in 
the perpendicular direction, corresponding to rotations 
around the tube axis, so that the amplitude is constant 
over the angular variable 0. For 0^0 the K-K' de- 
generacy is broken, and thus this angular distribution is 
expected. 

The eigenstates do not have a well-defined spin because 
of the SOI. However we can refer here to a "spin-up" 
state because the amplitude of the spin down component 
is several orders of magnitude smaller than that of the up 




FIG. 9: a) Density of states of a chiral (12,9) CNT (ry = 25°) 
with a length of 16 unit cells (about 42 nm) taking into ac- 
count curvature effects, the spin-orbit coupling and the Zee- 
man effect, b) Density of states of an armchair CNT (10,10) 
of 172 unit cells (about the same length of 42 nm) and with 
the same approximations as in (a), c) and d) are the dl/dV 
characteristics in units of Go corresponding respectively to a) 
and b). The value of the SOI 6 is set to 0.0028. Notice the 
absence of localization in the armchair case. 

This localization induced in the magnetic field is also ob- 
served in chiral nanotubes as mentioned before in Sec. III. 
We can see in Fig. 9a,c the case of a (12,9) chiral CNT, 
which has a chiral angle already quite close to the one 
of armchair nanotubes, and for which the localization 
is clearly visible. For armchair nanotubes on the other 
hand (Fig. 9b, d) this localization is not possible, as the 
boundary conditions do not allow for imaginary solutions 
of the K vector. 

The very good agreement between numerical and analyt- 
ical results is clearly seen from Fig. 10. It makes evident 
that the spin-flipping part of the Hamiltonian has a very 
small contribution so that we are justified to assign a 
specific spin to the transport features. In this way we 
can see that it is possible to have a control of the valley 
and spin degrees of freedom in CNTs due to the presence 
of SOI. This feature will thus be exploited in the next 
section. 
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FIG. 10: Comparison of numerical and analytical results for 
the projected density of states in different CNTs of approxi- 
mately the same length (about 42 nm). The dashed lines in 
the analytical calculation correspond to states in the K cone 
and the solid lines correspond to states. Spin up are rep- 
resented in red, and spin down in blue. All the lines shown 
correspond to the band with /^^ = 0, except for the states 
appearing in the zigzag CNT at the Fermi energy at = 0, 
which are the edge states corresponding to neighboring bands 
(the next allowed hi± values, hi± — i;^) and for which dotted 
lines are used, a) A (9,0) CNT as in Fig. 6d. b) The same 
(9,0) CNT but with a more realistic value of ^, as in the cal- 
culations shown in Fig. 9. c) A (12,9) CNT as in Fig. 9a. d) 
A (10,10) CNT as in Fig. 9b. 



VI. GENERATION OF SPIN-POLARIZED 
CURRENTS 




FIG. 11: a) Spin-resolved energy spectrum for the CNT 
(12,9) with a length of 100 unit cells and a SOI parameter 
b of 0.013. The colors correspond to those in Fig. 10. b) 
Corresponding differential conductance, where the spectrum 
lines have been symmetrized due to the considered symmet- 
rical voltage drop. 

The presence of localized states of definite spin, suggests 
the possibility of generating polarized spin currents, with 
in addition a definite isospin, by exploiting the localiza- 
tion phenomenon. For example, if a spin up current is 
desired, one should look for CNTs such that the local- 
ization onset of the K ^ states occurs already at neg- 




0.01 0.02 0.01 0.02 




FIG. 12: Spin-resolved energy spectrum for (12,9) CNTs of 
several lengths with a SOI characterized by the parameter 
8 = 0.0028. a) A hundred unit ceh long CNT, that is, about 
259 nm long, b) A CNT with 200 unit ceUs, and thus a length 
of about 518 nm. c) CNT with a length of about 1295 nm (500 
unit cells), d) Evolution with the magnetic field of the band 
edges for an infinite CNT. The colors correspond to those in 
Fig. 10. 



ative magnetic fields. This condition would lead to a 
region around the Fermi energy of highly spin-up polar- 
ized current for positive magnetic fields. Correspond- 
ingly spin-down polarized current characterizes negative 
fields. This possibility is demonstrated in Fig. 11a with 
^^oc(K|) = —0.001^0 7 calculated for clarity with a large 
value of 5. In a wide energy window around the Fermi 
energy all transport features are characterized by a spin 
up, with the only exception of the states which are 
governed by wave functions decaying already exponen- 
tially from the edge states as seen in Fig. lib. We can 
achieve this state also with a realistic ^ by tuning the 
other parameters influencing (lyioo as the length of the 
CNT. An example of a CNT also fulfilling these condi- 
tions is for instance a (38,37) CNT with a length of 100 
unit cells. Its chirality is very close to that of armchair 
CNTs, which minimizes the shift A/c^, and together with 
its long length allows for a localization onset of K | at 
negative fields. A greater length is of advantage because 
for larger CNTs the onset of the localization is moved 
towards smaller magnetic fields, and thus corresponds to 
a smaller Zeeman shift. We can see this effect in Fig. 12, 
where the spectra for (12,9) chiral CNTs with different 
lengths are shown. The features of the resonances evolve 
more and more towards the results expected for the in- 
finite CNTs (Fig. 12d), but showing the localization in- 
duced by the magnetic field at finite lengths. 

Another kind of polarization shows up in these systems as 
a result of the combination of the SOI and the localized 
states in the magnetic field. Without the SOI we have 
observed in Fig. 6c that the localizing states cross the 
Fermi energy in the presence of the Zeeman effect at the 
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FIG. 13: a) Evolution of the spectrum in the magnetic field 
and b) differential conductance as well as c) integrated con- 
ductance for a (12,9) CNT of 200 unit cells. The points 
around Vh — have been suppressed in order to avoid a di- 
vision by zero, d) Growth with the CNT length of the sepa- 
ration in magnetic flux of the conductance peaks around the 
Fermi energy. 



same field intensity for the K | state of the conduction 
band and the K \ state of the valence band. At this 
crossing point we have therefore a mixture of spin-up 
and spin-down states. The inclusion of SOI breaks the 
electron-hole symmetry of this picture and this crossing 
point is split into two, as the state lying closer to the 
Fermi energy crosses it at a lower magnetic field than the 
state. In a small bias window we thus obtain a spin- 
down polarized current swapping into a spin-up polarized 
current with increasing magnetic field. These two current 
peaks can be resolved in the magnetic flux as seen in 
Fig. 13. For negative magnetic fields both the spin and 
the isospin of these states will be reversed. Fig. 13 shows 
a CNT of 200 unit cells as the one described in Fig. 12b. 
The distance between these peaks of polarized current is 
set by the SOI parameter as well as the characteristics 
of the tube, such as length, radius and chirality. This 
separation of the peaks in the magnetic field will increase 
with the CNT length as seen in Fig. 13d. 



VII. CONCLUSIONS 

We have calculated the energy spectrum and transport 
properties of finite carbon nanotubes under a parallel 
magnetic field, including spin-orbit interactions. The 
spin-orbit coupling was treated following the approach 
of Ando,^^ but the spin quantization axis has been cho- 
sen parallel to the magnetic field to conveniently incor- 
porate its effects. The mixing of the tt and cr orbitals, 
due to the curvature of the nanotube, is automatically 
taken into account by this approach. 
We derived an analytical formula for the energy spectrum 
of these systems, which provides a very good matching to 



the numerical calculations in its limits of validity, i. e. for 
large enough radii to justify the approximations made. 
For radii of about 0.4 nm the agreement is already over 
90%. This analytical formula can be applied to finite-size 
systems if the proper boundary conditions are imposed. 
The very good results of this approach allow us to use 
it as a powerful tool to gain insight into the numerical 
results, as we can identify the composition of the states in 
the energy spectrum without the highly costly numerical 
calculation of the eigenstates. We can furthermore vary 
the length of the central system at no computational cost. 
The analysis of the transport properties {e.g. the dl /dV 
characteristics) shows a mechanism of localization in chi- 
ral and zigzag nanotubes in a parallel magnetic field. The 
magnetic field modifies the boundary conditions, causing 
hitherto extended states to localize near the ends of the 
tubes. This localization is gradual and initially the states 
involved are still conducting. But at a critical value of 
the magnetic field, which depends on the nanotube chi- 
rality and length, the localization is complete and the 
transport is then suppressed. 

Furthermore the SOI together with the Zeeman effect 
introduces a spin-splitting that breaks the electron-hole 
symmetry and allows the separation of polarized states. 
By exploiting this effect we achieve spin polarized cur- 
rents across the system. The magnetic field can thus be 
used as a tool for controlling the spin and isospin of the 
current through the CNT. 
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Appendix A: Modification of the hopping integrals 
by the SOI 

In the lowest order perturbation theory the state \zj 'fy) 
is given by 

a=x,y,z ^2p 
as^zty 

(Al) 

The kets | '\y ) and | ) are the two eigenstates of the 
Pauli's matrix ay, the normalized eigenspinors in the yj 
direction which coincides with the global Y. They can 
be expressed as functions of the eigenspinors in the zj 
direction, the original axis of spin quantization, as: 

iej/2 

lt.) = -^{|t,)+iU.)}, 
U.) = ^^{|t,)-iU.)}. 
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The action of Pauli's spin matrices dx and cr^ on these 
spinors produces a spin flip whereas (jy conserves the 
spin. Because of the orthogonality of the spin eigenstates 
and the symmetry of the orbital parts of the integrals, 
we can see that most terms will vanish, leaving only two 
contributions: 



+ 1% iy ) 



( -^j '^y \^so\^j 'ty ) 
iVj iy l^sol^j t^y ) 



(A3) 



where = (VF x p)^ aa and e^^ = e^p - e^^ > 

is the energy difference between the Pz and sp'^ orbit als, 
which are building respectively the tt and a molecular 
orbitals. It can already be predicted by symmetry that 
the orbital contributions of the integrals are equal, so 
that we can define 

e^j = -i{xj iy W = ^~'^'{yj U \Ko\^j t^y), 

(A4) 

in order to write the modified 7r-states as 

W - \zj ty ) + iS\xj ty ) + Se'^^ \yj iy ), (A5) 
and similarly for spin down 

\Zj iy) ^ \Zj iy ) - iS\Xj iy ) " Sc''^^ \y j ^y ). (A6) 

We have arrived at an expression for the perturbed or- 
bitals in terms of a global spin direction. The parameter 
S contains then the measure of the strength of the intra- 
atomic spin-orbit coupling. 

We proceed now to calculate the modified hopping pa- 
rameters between these states, which, after carrying out 
the spin part of the integrals and discarding higher order 
terms in (5, are given by 

{Zi^y \n\z,ty) = izmzj) 

+i5{{zi\n\xj)-{xmzj)}, 
{zi iy I H \zj iy) = ( ZilHlzj ) 

-is{{z,\n\xj)-ixi\n\zj)}, 

{Ziiy \n\zjiy) = 5e-''*{yi\n\zj) 

-Se-''^{zm\yj), 

{zUy \n\zjty) = 5e''^{zi\H\yj) 

-Se''^{yi\n\zj). (A7) 

The last two hopping integrals are the spin flipping terms, 
which the SOI make possible. They are a consequence 
of the spin mixing appearing in Eq. (A5) and Eq. (A6), 
where we see that the states including the SOI do not 
have a well-defined spin, but there is some admixture of 
the opposite spin species. 

In order to calculate the orbital part of the transfer inte- 
grals in Eq. (A7) we observe that the hopping {ai\l-L\aj) 
between the neighbor orbitals ) and |aj ), a = x, 
can be written in terms of its normal and tangential com- 
ponents in the graphene plane building the nanotube, as 



we have already seen in Sec. II B, for the case of \zj) 
orbitals. Explicitly, 



\Rji\ 



(A8) 



for generic Pa, orbitals. 
Thus we find 

( Zimx, ) = -( Ximzj ) = V;^ sin (0, - e,) 

-iVpp - VpD^si^ {Oi - ej) [cos {e, - e,) - 1] , 



izi\n\yj) = -iyi\'H\zj) 

= {v;^-v;^)-^[cos{e,-e,)-i], 

C 

i zi\n\zj) = v;^ cos (0,-0,) 

-{v;p-v;^)^[cos{o,-ej)-i]\ 



(A9) 



where Yji = Yj —Yi. This directly leads us to the final 
expression of the modified hopping parameters as given 
in Eq. (9). 



Appendix B: Approximations towards an effective 
Hamiltonian 



The transfer integral t^, which is a 2x2 matrix block as 
the spins are no longer degenerate, can be written as 



ZL 



^K^pp ^ppJor,^_ \ R J JIQ _^ 



-L „ ^9 



e-^^ 



(Bl) 



where we have defined 



2 , and df- and (i{ are 



the components of the vector in the basis 
of direction perpendicular and parallel to the tube axis, 
see Fig. Ic. The first row in the right-hand side of 
Eq. (Bl) contains only spin-independent elements, the 
second row groups spin-dependent terms whereas the last 
one is formed by a spin-flipping term. This last term de- 
pends on the local position of the interacting atoms in 
the tube through the angle 0. 

To obtain Eq. (Bl) from Eq. (9) we have taken the limit 
of large radii so that ac/i? is a quantity small enough for 
the following substitutions to be reasonable: 

• ^(R) ^ e As seen in Fig. lb this is the 

case for a small enough {Oi — Oj) angle. 
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• sin {Oi — Oj) ! 

cos {6i — Oj) 



2K2 



R ' 



We will further expand the Hamiltonian around the Dirac 
points, K and K', as we are interested in the physics 
taking place around the Fermi energy. The choice of K = 



-(-1,0) and 



47r 



-(1,0) allows us to write the 



two independent Fermi points in a compact form as rK 
with r = +1 for K and r = —1 for K^ For the expansion 
we substitute k — > rK + k being small enough to 



justify the expansion e 



ikd. 



rK d, 



(1 



In this way the matrix block 1-L\b of the Hamiltonian 



reads 



ace 



irSr] 



^pp 2 




where we have kept terms up to the second order in 
the variables ac/R, ^ and S and neglected higher order 
terms. 



^ P. Avouris, Z. Chen, and V. Perebeinos, Nature Nanotech- 

nology 2, 605 (2007), . . .0 . 1038/nnano . 2007 . 
^ S. Reich, C. Thomsen, and J. Maultzsch, Carbon nanotubes 

(Wiley- VCH, Berlin, 2004). 
^ K. Tsukagoshi, B. W. Alphenaar, and H. Ago, Nature 401, 

572(1999), doi: 10. 1038/44108. 
^ J. R. Kim, H. M. So, J.-J. Kim, and J. Kim, Phys. Rev. B 

66, 233401 (2002), doi : 10 . 1103/PhysRevB . 66 . 233401. 
^ S. Sahoo, T. Kontos, J. Furer, C. Hoffmann, M. Graber, 

A. Cottet, and C. Schonenberger, Nature Physics 1, 99 
(2005), doi:10.1038/nphysl49. 

^ A. Cottet, T. Kontos, S. Sahoo, H. T. Man, M.-S. Choi, W. 

Belzig, C. Bruder, A. F. Morpurgo, and C. Schonenberger, 

Semiconductor Science and Technology 21, S78 (2006), 

doi:10.1088/0268-1242/21/ll/Sll. 
^ S. Koller, L. Mayrhofer, and M. Grifoni, New Journal of 

Physics 9, 348 (2007), doi : 10 . 1088/1367-2630/9/9/348. 
^ C. Schenke, S. Koller, L. Mayrhofer, and M. 

Grifoni, Phys. Rev. B 80, 035412 (2009), 

doi : 10 . 1 103/PhysRevB . 80 . 035412. 
^ F. Kuemmeth, S. Hani, D. Ralph, and P. McEuen, Nature 

452, 448 (2008), doi : 10 . 1038/nature06822. 
^° T. Ando, J. Phys. Soc. Japan 69, 1757 (2000), 

doi : 10 . 1 143/ JPS J . 69 . 1757. 

L. Chico, M. P. Lopez-Sancho, and M. C. 
Mufioz, Phys. Rev. Lett. 93, 176402 (2004), 
doi : 10 . 1 103/PhysRevLett . 93 . 176402. 
D. Huertas- Hernando, F. Guinea, and A. 
Brataas, Phys. Rev. B 74, 155426 (2006), 
doi : 10 . 1 103/PhysRevB . 74 . 155426. 

H. O. H. Churchill, F. Kuemmeth, J. W. Har- 
low, A. J. Bestwick, E. I. Rashba, K. Flensberg, 
C. H. Stwertka, T. Taychatanapat, S. K. Watson, 
C. M. Marcus, Phys. Rev. Lett. 102, 166802 (2009), 
doi : 10 . 1 103/PhysRevLett . 102 . 166802. 
S. H. Jhang, M. Marganska, Y. Skourski, D. Preusche, 

B. Witkamp, M. Grifoni, H. van der Zant, J. Wosnitza 
and C. Strunk, Phys. Rev. B 82, 041404(R) (2010), 
doi : 10 . 1 103/PhysRevB . 82 . 041404. 



D. V. Bulaev, B. Trauzettel, and D. Loss, Phys. Rev. B 
77, 235301 (2008), doi. xv.. xx03/PhysRevi^. i i 
T. S. Jespersen, K. Grove-Rasmussen, J. Paaske, K. 
Muraki, T. Fujisawa, J. Nygard, and K. Flensberg, 
Nature Physics, advanced online publication (2011), 
doi:10.1038/nphysl88 . 

H. Ajiki and T. Ando, J. Phys. Soc. Japan 62, 1255 (1993), 
doi : 10 . 1143/JPSJ . 62 . 1255. 

W. Tian and S. Datta, Phys. Rev. B 49, 5097 (1994), 
doi : 10 . 1103/PhysRevB . 49 . 5097. 

A. Bachtold, C. Strunk, J.-P. Salvetat, J.-M. Bonard, L. 
Forro, T. Nussbaumer, and Christian Schonenberger, Na- 
ture 397, 673 (1999), doi : 10 . 1038/17755. 

A. Fujiwara, K. Tomiyama, H. Suematsu, M. Yu- 
mura, and K. Uchida, Phys. Rev. B 60, 13492 (1999), 
doi : 10 . 1103/PhysRevB . 60 . 13492. 

S. Zaric, G. N. Ostojic, J. Kono, J. Shaver, V. C. Moore, 
M. S. Strano, R. H. Hauge, R. E. Smalley, and X. Wei, Sci- 
ence 304, 1129 (2004), doi : 10 . 1126/science . 1096524. 
U. C. Coskun, T.-C. Wei, S. Vishveshwara, P. M. 
Goldbart, and A. Bezryadin, Science 304, 1132 (2004), 
doi : 10 . 1126/science . 1096647. 

C. Strunk, B. Stojetz, and S. Roche, Semicon- 
ductor Science and Technology 21, S38 (2006), 
doi : 10 . 1088/0268-1242/21/11/S06. 

B. Stojetz, S. Roche, C. Miko, F. Triozon, L. Forro, 
and C. Strunk, New Journal of Physics 9, 56 (2007), 
.01 : 10 . l088/1367-26dO/y/3/Ov . 

G. Fedorov, P. Barbara, D. Smirnov, D. Jimenez, 
and S. Roche, Appl. Phys. Lett. 96, 132101 (2010), 
doi: 10. 1063/1.336021^. 

Y. Zhang, Q. Liang, and J. Dong, Physics Letters A 372, 
6996 (1998), qqx ; lu . 1016/j .physleta. 2008 . 10 . 007. 
S. Krompiewski, R. Gutierrez, and G. Cu- 
niberti, Phys. Rev. B 69, 155423 (2004), 
doi : 10 . 1103/PhysRevB . 69 . 155423. 
2^ N. Nemec and G. Cuniberti, Phys. Rev. B 74, 165411 
(2006), Hoi -10 1103/PhvsRevB.74.l65^ . 
S. Krompiewski and G. Cuniberti, Journal of Mag- 



14 



netism and Magnetic Materials 310, 2439 (2007), 
doi: 10. 1016/1 . immm. 2006 . 10 . 90 . 

W. Izumida, K. Sato, and R. Saito, J. Phys. Soc. Japan 
78, 074707 (2009), uo±: . ^j.^o/ Ji^bj . iq .ui^iu i . 
B. Wunsch, Phys. Rev. B 79, 235408 (2009), 
doi : 10 . 1 103/PhysRevB . 79 . 23540P. 

A. Secchi and M. Rontani, Phys. Rev. B 80, 041401 (R) 

(2009) , doi : 10 . 1 103/PhysRevB . 80 . 041404 . 

J. Zhou, Q. Liang, and J. Dong, Phys. Rev. B 79, 195427 

(2009) , doi : 10 . 1103/PhysRevB . 79 . 195427. 

J.-S. Jeong and H.-W. Lee, Phys. Rev. B 80, 075409 

(2009) , doi : 10 . 1 103/PhysRevB . 80 . 075409. 

T.-F. Fang, W. Zuo, and H.-G. Luo, Phys. Rev. Lett. 101, 

246805 (2008), ioi : 10 . 1103/PhysRevLett . 101 . 246805. 

D. E. Logan and M. R. Galpin, J. Chem. Phys. 130, 224503 

(2009), ..v,...^^._. 3148035. 

S. Weiss, E. L Rashba, F. Kuemmeth, H. O. H. Churchih, 
and K. Flensberg, Phys. Rev. B 82, 165427 (2010), 
doi : 10 . 1 103/PhysRevB . 82 . 16542' . 

M. Margahska, M. del Valle, S. H. Jhang, C. Strunk, and 
M. Grifoni, arXiv:1101.2088vl. 

Introducing molecular electronics, edited by G. Cuniberti, 
G. Fagas, and K. Richter (Springer, Berlin and Heidelberg, 
2005), doi:10.1007/bl01525. 

K. Sasaki, S. Murakami, R. Saito, and Y. 
Kawazoe, Phys. Rev. B 71, 195401 (2005), 
doi : 10 . 1 103/PhysRevB . 71 . 195401. 

We keep the conventional notation for the chirality of a 
CNT , so that with our choice of lattice vectors a nanotube 
of chirality (n,m) has a chiral vector C = (n + m)ai +ma2. 
R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physi- 
cal properties of carbon nanotubes (Imperial College Press, 
London, 1998). 

J. Luttinger, Phys. Rev. 84, 814 (1951), 
doi: 10. 1103/PhysRev.84.81^ . 

A is canceled out leaving us with the original Hamiltonian 
for zero magnetic field and no extra terms if assuming a 
spatially slow varying magnetic field and localized atomic 
functions. 

R. Peierls, Zeitschrift fiir Physik 80, 763 (1933), 
doi : 10 . 1007/BF01342591. 



H. Min, J. E. HiU, N. A. Sinitsyn, B. R. Sahu, L. Kleinman, 
and A. H. MacDonald, Phys. Rev. B 74, 165310 (2006), 
doi : 10 . 1103/PhysRevB . 74 . 165310. 

R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, Berlin, 
2003). 

P. R. Wallace, Phys. Rev. 71, 622 (1947), 
doi : 10 . 1103/PhysRev . 71 . 622. 

A. R. Akhmerov and C. W. J. Beenaker, Phys. Rev. B 77, 
085423 (2008) , doi : 10 . 1 103/Phy sRevB . 77 . 085423. 
L. Mayrhofer and M. Grifoni, Eur. Phys. J. B 63, 43 
(2008), doi : 10 . 1140/epjb/e2008-00204-0. 

We impose the wave functions ip± to vanish at the first row 
of missing atoms at the ends so that L = |T| {Nuc + 1/2), 
with T being the translational vector which has a length 
of v^ac for an armchair CNT, and Nuc being the number 
of unit cells in the tube. 

G. Cuniberti, F. Grossmann, and R. Gutierrez, 
Adv. Solid State Phys. 42, 133 (2002), 
doi : 10 . 1007/3-540-45618-Xll. 

D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981), 
doi : 10 . 1103/PhysRevB . 23 . 6851. 

D. Tomanek and S. G. Louie, Phys. Rev. B 37, 8327 (1988), 
ioi : 10 . 1103/PhysRevB . 37 . 8327 . 

E. D. Minot, Y. Yaish, V. Sazonova, and P. L. McEuen, 
Nature 428, 536 (2004), doi : 10 . 1038/nature02425. 

M. Fujita, K. Wakabayashi, K. Nakada, and K. 
Kusakabe, J. Phys. Soc. Japan 65, 1920 (1996), 
doi : 10 . 1143/ JPS J . 65 . 1920. 

K. Nakada, M. Fujita, G. Dresselhaus, and M. S. 
Dresselhaus, Phys. Rev. B 54, 17954 (1996), 
doi : 10 . 1103/PhysRevB . 54 . 17954. 

K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999), doi : 10 . 1103/PhysRevB . 59 . 8271. 
P. Kim, T. W. Odom, J.-L. Huang, and C. M. 
Lieber, Phys. Rev. Lett. 82, 1225 (1999), 
ioi : 10 . 1103/PhysRevLett . 82 . 1225. 

S. Okada and A. Oshiyama, J. Phys. Soc. Japan 72, 1510 
(2003) , doi : 10 . 1143/ JPS J . 72 . 1510. 



